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Abstract 
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— The  first  and  second  moments,  i.e.  the  coherent  field  and  the  two-point 
two-time  correlation  function,  are  calculated  for  the  acoustic  fields 


scattered  frcm  various  rough  surfaces.  For  each  surface  -they  yield  the 
reflection,  absorption  and  differential  scattering  coefficients,  as  well 
as  an  equivalent  boundary  condition  for  the  coherent  field.  Renormalized 
coefficients  are  constructed  to  eliminate  divergences  at  grazing  incidence. 
The  results  are  specialized  to  surfaces  which  are  statistically  homogeneous 
in  both  space  and  time,  to  surfaces  which  are  not  moving,  to  surfaces 
which  are  simply  or  multiply  periodic,  and  to  surfaces  consisting  of 
randomly  placed  bosses  on  a  smooth  surface.  The  surfaces  considered  are 
slightly  rough,  moving,  soft  or  hard  boundaries,  and  flat  surfaces  with  random 
admittances  or  impedances.  The  analysis  is  based  on  the  regular  perturbation 
method  or  Bom  approximation.  Comparisons  with  previous  results  are  made. 
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1.  Introduction. 


When  a  wave  is  incident  upon  a  rough  surface,  it  produces  reflected, 
transmitted,  and  scattered  waves.  The  energy  flux  in  each  wave  can  be 
determined  by  the  transport  theory  of  radiation.  This  theory  uses  geo¬ 
metrical  acoustics  to  describe  energy  propagation,  together  with  phe¬ 
nomenological  coefficients  to  describe  reflection,  transmission  or 
absorption,  and  differential  scattering  at  the  surface.  However  it  does 
not  determine  these  coefficients  in  terms  of  the  properties  of  the  sur¬ 
face  and  it  does  not  describe  the  fluctuation  and  correlation  properties 
of  the  waves. 

Several  methods  of  analysis,  based  upon  the  wave  equation,  have  been 
devised  to  overcome  these  two  shortcomings.  We  shall  employ  one  of  them, 
the  regular  perturbation  method  or  Born  approximation,  to  calculate  the 
first  and  second  moments  of  the  field.  These  moments  represent  the 
coherent  field  and  the  two-point  two-time  correlation  function  of  the 
field,  so  they  yield  information  about  cield  fluctuation  and  corre¬ 
lations.  From  them  we  shall  obtain  expressions  for  the  various  coeffi¬ 
cients  which  occur  in  transport  theory.  In  addition,  we  shall  show  that 
the  two-point  two-time  correlation  function  of  the  field  can  be  expressed 
in  terms  of  these  same  coefficients. 


Since  the  Born  approximation  has  been  used  to  analyze  reflection  by 
rough  surfaces,  the  present  work  may  be  viewed  as  a  completion  of  that 
analysis.  It  treats  in  a  unified  way  four  different  cases:  a  nearly 
flat  Soft  rough  surface,  a  nearly  flat  hard  rough  surface,  a  flat  surface 
with  a  small  random  impedance,  and  a  flat  surface  with 


a  small  randan  admittance.  The  rough  surface  is  allowed  to  have  a 
small  random  velocity,  and  the  random  admittance  and  im¬ 
pedance  are  allowed  to  vary  in  time.  The  first  and  second  moments  of  the 
field  are  obtained  in  all  these  cases,  which  has  not  been  done  before. 
Furthermore  the  coefficients  of  reflection,  etc.,  are  modified  or  re¬ 
normalized  to  eliminate  divergences  which  occur  at  grazing  incidence. 

From  the  general  results,  more  explicit  ones  are  obtained  for  a 
number  of  special  kinds  of  surface  roughness.  First  surfaces  which  are 
statistically  homogeneous  or  stationary  in  space  and  time  are  considered. 
Then  surfaces  composed  of  bosses  randomly  located  on  a  flat  surface  are 
treated.  These  are  the  surfaces  which  have  been  studied  extensively  by 
Twersky  [  1  ]  using  a  very  different  method  of  analysis,  so  we  compare 
our  results  with  his.  In  this  way  we  are  able  to  relate  the  results  of 
the  perturbation  method  to  those  of  Twersky's  self-consistent  field  method. 

A  survey  of  the  theory  of  scattering  from  rough  surfaces  together  with 
references,  is  contained  in  pages  54-73  of  the  review  article  by  DeSanto  [2], 


2.  Formulation  and  solution  in  the  nonrandom  case. 


Let  *Kx,y, z,t)  and  <J>(x,y,z,t)  be  the  incident  and  scattered  acoustic 
velocity  potentials  in  a  uniform  medium  of  density  p  and  sound  speed  c. 

We  assume  that  i|i  is  a  given  incoming  solution  of  the  wave  equation  near  a 
bounding  surface,  while  everywhere  above  this  surface  4>  is  an  outgoing 
solution  of  the  wave  equation: 

(A  -  c'232)<f>  =  0  .  (1) 

We  shall  consider  four  different  boundary  value  problems  for  (1)  and  de¬ 
note  them  a,  b,  c,  d.  The  first  two  deal  with  a  nearly  flat,  moving  sur¬ 
face  z  =  eh(x,y,t).  Here  e  is  a  small  parameter  and  h  is  a  given  func¬ 
tion. 

Case  a  represents  a  soft  or  pressure  release  surface  on  which  the 
pressure  vanishes, 

p3t0l>  +  <J>)  *  0  on  z  »  eh  .  (2af) 

Case  b  represents  a  hard  or  rigid  surface  on  which  the  normal  component 
of  acoustic  velocity  equals  the  normal  velocity  of  the  surface, 


Oz  -  ehx3x  -  ehy3y)0J>  ♦  ♦)  =  eht 

on  z  *  eh  . 


(2b*) 


The  right  side  of  (2b*)  is  an  inhomogeneous  or  source  term  for  the  radia¬ 
tion  produced  by  the  moving  surface.  We  shall  not  consider  this  radia¬ 
tion,  so  we  shall  investigate  the  homogeneous  form  of  (2b*),  which 
corresponds  to  a  vanishing  normal  component  of  acoustic  velocity  on  the 
surface . 


In  addition  to  non-flat  surfaces,  we  shall  consider  the  flat  surface 


z  =  0.  Case  c  involves  the  admittance  boundary  condition 

(eh3t  +  c3  )(ip  +  <p)  =  0  on  z  =  0  .  (2c) 

t  z 

In  case  d  we  have  the  impedance  boundary  condition 

(3  +ehc3  )0I>  +  <fr)  *  0  on  z  =  0  .  (2d) 

t  z 

The  specific  admittance  in  (2c)  and  the  specific  impedance  in  (2d)  are 
both  equal  to  eh. 


To  solve  these  four  problems  we  first  expand  <t>  +  ip  in  (2a'),  and  in 
the  homogeneous  form  of  (2b'),  in  a  Taylor  series  in  z  about  z  =  0,  and 
set  z  =  eh  to  obtain 


[1  *  eh3  + 
1  z 

P,  *  e(h^ 


2  2 

eV  ~2 

2  z 


-  h  3 
x  x 


+  0(e3)]3  ($  +  i|»)  *  0  on  z  *  0  , 

v 

2 

-  h  3  )  +  e2(V  33  -  hh  3  3  -  hh  3  3) 

y  y  2  z  xxz  y  y  z 

*  0(e3)  ]  (<J>  +  ip)  -  0  on  z  *  0  . 


(2a) 


(2b) 


Next  we  introduce  the  vectors  p  =  (x,y,t)  and  s  =  (ot,3,a>)  and  the 
Fourier  transform  pair 


F(s)  = 

(2ir)-3  j  f(p)e"ip*sdp  , 

(3) 

f(p)  = 

|  F(s)eip*Sds  . 

(4) 

We  shall  use  the  corresponding  capital  letter  to  denote  the  transform  of 
a  function  represented  by  a  small  letter.  When  f(p)  is  real  we  have 
F(-s)  ■  F(s)  where  the  overbar  denotes  the  complex  conjugate. 


Upon  applying  the  transform  (3)  to  (1)  we  find  that  $(s,z),  the 


The  solution  of  (5)  which  is  bounded  at  z  *  when  k  is  imaginary, 

or  which  is  outgoing  when  k  is  real,  is  of  the  form 

*(s,z)  *  As(s)e~ik(s)*  .  (8] 

We  assume  that  the  transform  T(s,z)  of  the  incident  field  has  the  form 

T(s,z)  -  Ai(s)eik(s)z  .  (9) 

From  (8)  and  (9)  we  see  that  the  propagating  modes,  those  for  which 
2  2  2 

(<u/c)  >  a  +  8  ,  are  incoming  in  f  and  outgoing  in  $.  Similarly  the 

evanescent  modes  in  $  decay  with  increasing  distance  z  from  the  boundary, 
while  those  in  f  grow  with  z.  Therefore  the  representation  (9)  of  V  can¬ 
not  hold  for  z  arbitrarily  large  if  it  contains  evanescent  modes. 

The  amplitude  A^(s)  in  (9)  is  supposed  to  be  known  while  the  ampli¬ 
tude  As(s)  in  (8)  is  to  be  found.  To  find  it  we  apply  the  transform  (3) 


to  the  boundary  conditions  (2a,b,c,d)  and  then  use  (8)  and  (9)  in  the  re¬ 
sulting  equations.  In  each  case  this  leads  to  an  integral  equation  for 


Ag(s).  To  solve  it  we  write  Ag(s)  in  the  form 

As(s)  =  A<0)(s)  +  eA^Cs)  ♦  e  V2)  (s)  +  0(e3)  .  (10) 


Then  we  use  (10)  in  the  integral  equation,  equate  coefficients  of 
corresponding  powers  of  e,  and  obtain  equations  for  the  determination  of 
Af>,  A^  and  Ag2^.  Next  we  solve  these  equations,  as  is  done  in 
Appendix  A  for  the  soft  boundary  (case  a) . 

We  can  write  the  result  (10)  for  Ag(s)  in  all  four  cases  in  the  form 
of  a  suitable  linear  operator  L  acting  on  A^.  Explicitly  we  find 


As(s)  =  [LA.](s)  =  +  (Ai(s) 

+  2e  |  HCs-s^m^s'M^s^ds' 

♦  2e2  J  iKs-sOn^s.s') 


(11) 


H(s,-s,,)ra(s^s,,)Ai(s,,)dsMds, 


0(e3)} 


In  (11),  H  is  the  transform  of  h.  The  function  m 

and  the  choice  of  sign  are  listed  in  Table  I  for  the  four  cases.  This 
result  (11)  gives  the  Fourier  transform  of  the  amplitude  of  the  reflected 
wave  for  the  deterministic  problem. 


3.  Reflection  coefficients  in  the  random  case. 


We  shall  now  use  the  preceding  results  to  treat  the  case  in  which 
the  function  h(x,y,t),  describing  the  admittance,  impedance,  or  the  sur¬ 
face,  is  a  sample  function  of  a  random  process.  This  means  that  h 
depends  upon  a  parameter  which  is  a  random  variable  having  an  associated 
probability  distribution.  Consequently,  functionals  of  h,  such  as  H,  cj> , 

$ ,  and  Ag,  are  random  and  depend  upon  the  same  random  variable.  We  per¬ 
mit  the  incident  field  ip,  and  hence  f  and  to  be  either  deterministic 
or  random.  However,  for  the  random  case,  we  restrict  our  attention  to  \p 
and  h,  and  therefore  A^  and  H,  which  are  statistically  independent. 
Without  loss  of  generality,  we  proceed  as  if  ij>  and  A^  were  random,  since 
the  deterministic  case  is  the  zero  fluctuation  limit  of  the  random  case. 

The  expected  value  or  statistical  mean  value  of  a  random  field  is 
called  the  coherent  field.  We  use  angular  brackets  to  denote  the  mean  or 
ensemble  average.  For  example,  <$>  is  the  coherent  scattered  field.  Its 
transform  «&>,  obtained  from  (8)  and  (11)  by  averaging,  is 


<$(s,z)>  =  <As(s)>e"xk(s)z 

=  [<L><Ai>](s)e"lk(s)z  . 


(12) 


The  average  operator  <L>,  obtained  by  averaging  (11),  transforms  a  func¬ 
tion  F  according  to  the  rule 


[<L>F]  (s)  =  +  {F(s) 


+2 e  J  ^(s-s'^m^s'JFfs^ds’ 

( 

+  2e2  |J  <H(s-s,)H(s,-sM)>m(s,s,)m(s',s,')F(s,,)ds,ds" 


+  0(e)}  . 


We  note  that  only  second  order  moments  of  H  and  h  are  required  to  compute 

3 

the  coherent  field  when  we  neglect  0(e  ). 

The  average  of  the  scattered  amplitude  simplifies  significantly  if 
the  process  h  is  statistically  second  order  stationary  in  space  and  time. 


This  means  that  <h>  is  a  constant,  hQ,  and  the  two-point  correlation  r  of 
depends  only  on  the  difference  in  arguments.  Therefore,  we  have 


<h(p)>  =»  hQ  ,  (14) 

<h(p  +  p')h(p')>  =  r(p)+  h2  .  (15) 

These  equations  imply  that  the  first  two  moments  of  the  transform  H  are: 

<H(s)>  =  h0<5(s)  ,  (16) 

<H(s)H(s')>  *  R(s)6 (s  +  s*)+  h26(s)6(s')  . 


(17) 


10 


Here  R(s),  the  spectral  power  density  of  h,  is  the  Fourier  transform  of 
the  auto-correlation  function  r,  and  5  is  the  Dirac  delta  function.  The 
function  R(s)  is  a  positive,  even,  real-valued  function  of  s. 
Substituting  (16)  and  (17)  into  (13)  yields: 

<L>  =  C(s) I  +  0(e3)  .  (18) 


Here  I  is  the  identity  and  C(s)  is  given  by 
C(s)  =  +jl+  2eh0m(s,s)  +  2£2h(2in2(s,s) 

+  2e2  |  R(s-s,)m(s,s,)m(s’ ,s)ds' +  0(e3) 

From  (12)  and  (18),  we  see  that  C(s)  is  the  ratio  of  the  coherent 
scattered  amplitude  to  the  coherent  incident  amplitude 


(19) 


<As(s)>  =  C(s)<A.(s)>  +  0(e3)  . 


(20) 


Thus  C(s)  is  the  reflection  coefficient  of  the  surface.  It  is  given  in 
Table  II  for  the  four  cases,  based  upon  (19)  and  Table  I.  If  the  incident 
wave  is  a  plane-wave,  (20)  holds  for  the  actual  wave  amplitudes,  as  well  as 
for  the  transformed  amplitudes.  Hence,  C(s)  is  also  called  the  plane  wave 
reflection  coefficient. 


11 


4.  Second  moments  of  the  acoustic  potential. 

We  now  use  the  previous  results  to  evaluate  y>  the  second  moment  of 
the  total  acoustic  potential.  This  two-point,  two-time  autocorrelation 
of  +  <p  is  defined  by 

Y(x,X)  =  <[V>  +  <J>](X  +  x/2)[ip  *  <fr](X  -  x/2)>  .  (21) 

Here  X  ■  (X,Y,Z,T)  is  the  midpoint  of  the  two  points  at  which  y  is 
evaluated  and  x  =  (x,y,z,t)  is  the  difference  between  the  two  points. 

The  transform  of  y  with  respect  to  (x,y,t)  with  transform  variable 
s  ■  (a,S,w),  and  then  with  respect  to  (X,Y,T)  with  transform  variable 
q  =  (a,b,w)  is  denoted  by  T.  We  have  from  (8)  and  (9) 

r(s,q;z,Z)  =  <[Y  +  *](q/2  +  S,Z  +  z/2)  [?  +  *] (q/2  -  s,  Z  -  z/2)> 

*  Bii(s,q)expi{+  [k+  +  k_]Z  +  [k+  -  k_]z/2> 

+  B.  (s,q)expi{+  [k  -  k  ]Z  +  [k  +  k  ]z/2) 

is  +  "  +  (22) 

+  Bgi(s,q)expi{-  [k+  -  k_]Z  -  [k+  +  k_]z/2> 

+  Bss(s,q)expi{-  [k+  +  k_]Z  -  [k+  -  k_]z/2>  . 

Here  k+  =  k(q/2  +_  s) . 

In  (22)  the  airplitude  correlations  B^^,  etc.  are  defined  in  terms  of 
the  scattered  and  incident  amplitudes  as  follows: 

Bii(s,q)  =  <A.(q/2  +  s)A.(q/2  -  s)>  , 

Bsi(s,q)  =  <Ag(q/2  +  s)A.(q/2  -  s)>  =  <[<L>Ai](q/2  +  s)Ai(q/2  -  s)> 

=  Big(-s,q)  ,  (23) 

B  (s,q)  =  <A  (q/2  +  s)A  (q/2  -  s)>  =  <[<L>A  ]  (q/2  +  s)  [<L>A.  ]  (q/2  -  s)> 

SS  S  S  l  ^ 

♦  <[L»A.](q/2  +  s)[L'A.](q/2  -  s)>  . 


The  average  <L>  is  given  by  (13)  and  the  fluctuating  part  of  L,  denoted 
L'  =  L  -  <L>,  is 


► 

[L'F]  (s)  =  +_2e  [H(s-s')  -  ^(s-s'^lnKSjSMFCs^ds1 

4 


*  0(e2) 


Only  the  order  e  term  of  L*  enters  the  computation  of  B  ,  provided  we 
3 

neglect  0(e)  terms. 

The  expression  for  the  second  moment  simplifies  when  h  is  second 
order  stationary.  Then  (23)  yields 


Bsi(s,q)  =  C(q/2  +  sjB^s.q)  +  0(e3)  ,  (25) 

Bss(s,q)  =  C(q/2  ♦  s)C(q/2  -  s)B..(s,q) 

2 

+  4e  R(s-s,)m(q/2+s,q/2+s,)m(q/2-s,q/2~s,)Bii(s '  ,q)ds' 

+  °(e3)  •  (26) 


Further  simplification  results  if  the  incident  field  ^(x,y,z,t)  is 
statistically  second  order  stationary  in  the  variables  x,  y,  and  t.  For 
this  case,  we  have: 

Bu(s,q)  =  I.(s)6(q)  ,  (27) 

where  I^s)  is  the  intensity  of  the  incident  field.  With  (27),  equations 
(25)  and  (26)  become 

Bgi(s,q)  =  C(s)Ii(s)6(q)  +  0(e3)  ,  (28) 

Bss(s>q)  =  +  °Ce3)  • 


(29) 


Here  Is(s),  the  intensity  of  the  outgoing  field,  is  defined  by 

Is(s)  =  |C(s)]2Ii(s)  +  k  2(s)  J  aCs.s^I^s'Jds'  +  0(e3)  .  (30) 

In  (30),  we  have  introduced  the  important  quantity  a(s,s'),  the  differen¬ 
tial  scattering  coefficient: 

a(s,s')  =  4e2k2(s)R(s-s') |m(s,s')  |2  .  (31) 

The  scattering  coefficients  for  admittance,  impedance,  soft,  and  hard 
boundary  conditions  appear  in  Table  III.  Since  R  ^  0,  we  note  that  ex  ^  0. 

The  physical  interpretation  of  Ig(s)  is  evident  from  its  definition 
(30  )•  It  consists  of  the  incident  intensity  function  I^(s)  multiplied  by 
the  absolute  square  of  the  reflection  coefficient  for  wavevector  s,  plus 
the  integral  over  the  contributions  from  the  incident  waves  with  wave- 
vector  s'  which  are  scattered  from  s'  to  s.  Thus  it  represents  the  sum  of 
the  reflected  and  scattered  intensities  with  wavevector  s. 

By  using  the  results  for  stationary  processes  in  (  22 ) ,  and  again 
making  use  of  (  7),  we  obtain  the  following  form  for  T: 

r(s,q;z,Z)  =  5(q)[I.(s)eikCs:)z  +  C(s)I.(s)e2ik(s)Z 

+  C(s)I.(s)e"2ikCs)Z  *  Is(s)e'ik(s)z]  +  0(e3)  , 

2  2  2 

(oj/c)  >_  a  +  , 

(32) 

=  <S(q)[I.(s)e2ik(s)Z  +  C(s)I.(s)eikCs)z 

+  C(s)I  (s)e"ik(s)z  +  I  (s)e'2ik(s)Z]  ♦  0(e3)  , 

2  2  2 

(co/c )  <  a  ♦  B  . 


The  result  (32)  gives  the  Fourier  transform  of  the  two-point  two- 
time  correlation  function  of  the  total  potential  ip  +  <j>.  It  applies  when 
the  function  h  and  the  incident  field  \p  are  both  statistically  homogeneous 
with  respect  to  x,  y  and  t.  In  particular  it  holds  for  a  time-harmonic 
plane  wave  with  a  random  phase  which  is  uniformly  distributed  in  the  in¬ 
terval  (0,2ir).  It  also  holds  for  a  non-random  time  harmonic  plane  wave 
provided  that  the  averages  denoted  by  angular  brackets  are  interpreted  to 
mean  stochastic  averaging  and  time  averaging  with  respect  to  T.  From 
(32)  we  see  that  the  total  potential  is  second  order  stationary  with 
respect  to  the  midpoint  coordinates  X,  Y  and  T,  but  not  with  respect  to  Z. 

The  corresponding  result  for  the  scattered  field  alone  is  given  by 
the  term  in  (32)  proportional  to  Is>  Denoting  it  we  have 

T  (s,q;z,Z)  =  I  (s)6(q)e'lk(-s)z  +  0(e3)  ,  (to/c)2  >  a2  +  B2  , 

ss  S  “  (33 

=  Is(s)S(q)e"2ik(s)Z  ♦  0(e3)  ,  (w/c)2  <_  a2  +  B2  . 

This  is  the  result  for  the  Fourier  transform  of  the  two-point  two-time 
correlation  of  the  scattered  field  in  the  stationary  case.  We  see  from 
the  expression  (30)  for  Is(s)  that  in  addition  to  the  incident  intensity 
I. (s),  T  _  involves  only  the  magnitude  of  the  reflection  coefficient, 
|c(s)|,  and  the  differential  scattering  cross-section  <j(s, s' )  of  the 


surface. 


flux  and  energy  densit 


The  energy  flux  vector  j (X)  and  the  energy  density  e(X)  are  defined 
by 


j  (X)  *  -p3T0J>  +  4>)*vx0|>  +  <t>)  ,  (34) 

e(X)  =  |  {c‘2[3t(4)  +  4>)]2  +  +  4>)]2}  .  (35) 

From  the  wave  equation  (1)  satisfied  by  ip  +  <p  in  a  source-free  region,  it 
follows  that  j  and  e  are  related  by  the  energy  conservation  equation 

3t?  +  7X*j  =  0  *  (36) 

When  ip  +  4>  is  random,  the  average  flux  <j>  and  the  average  energy 
density  <e>  are  defined  by  averaging  both  sides  of  (34)  and  (35) .  These 
average  quantities  then  satisfy  the  averaged  energy  conservation  equation 
(36).  They  can  be  expressed  in  terms  of  the  two-point  two-time  correla¬ 
tion  function  y(x,X) .  This  is  shown  for  <j>  by  the  following  chain  of 
equalities: 

<j(X)>  =  -p{<3T[4>  +  4>](x)»Vx[i|>  +  4>](x)>>x=x 
=  -pOTVxY(X  -  x,j[X  ♦  x]))x=x 

*»t<St7x-  J  3T7XWx-X)1x=0  •  l37) 


The  last  equality  of  (37)  is  a  consequence  of  the  evenness  of  y  in  the 
separation  vector  x.  A  similar  analysis  yields  the  energy  density  at  the 
point  X: 
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<e(X)>  =  -  |  p([c_2(32  -  is2)  +  (Ax  -  i  Ax)]y(x,X)}x=0  .  (38) 

When  h  and  <p  are  statistically  stationary  in  x,  y  and  t,  then  y  is 
independent  of  X,  Y  and  T,  as  we  have  just  shown  in  the  previous  section. 
Therefore  from  (37)  and  (38)  it  follows  that  <j>  and  <e>  are  functions 
only  of  Z.  Then  the  averaged  form  of  (36)  reduces  to  3  <  j  >  =  0,  which 
shows  that  <  j  >,  the  z  component  of  <j>,  is  constant.  By  using  (37) 
for  <j>  and  evaluating  the  right  side  at  X  =  0,  we  can  write  this  result 
as 


<jzCX)> 


p3t3zy(0,0) 

* 

p  iai3zr(s,q,0,0)dsdq  . 


(39) 


The  second  equality  follows  from  the  definition  of  the  transform  T  of  y 
at  X  =  0. 

We  now  use  (32)  for  T  in  (39)  and  then  use  (30)  to  express  1^ .  Then 
we  can  write  <jz>  in  the  form 


<j2(X)> 


a(s)pwk(s) I^fsjds  . 


(40) 


Here  we  have  introduced  the  absorption  coefficient  a(s)  defined  by 


a(s)  =  1  -  | C(s) 


^k(s)k(s' ) * s)ds'  +  0(e  )  * 

2  2  2 
(co/cr  >  a  ♦  8  , 


(41) 


=  C*(s)  -  C(s)  -7^py]a(s’,s)ds'  +  0(e) 


2  2  2 
(ai/c)  <  a  ♦  S 


The  integrals  are  over  the  propagating  modes,  (co'/c)  >ot  +  S  ,  as  is 
indicated  by  the  symbol  >.  For  propagating  modes,  (41)  shows  that  a(s) 
is  that  fraction  of  the  incident  energy  flux  which  is  neither  reflected 
nor  scattered. 


6.  Renormalization. 

Some  of  the  main  results  we  have  derived  contain  singularities  or 
divergences  which  are  physically  unreasonable  and  mathematically  inappro¬ 
priate.  Thus  for  example,  the  reflection  coefficient  C(s),  the  scattered 
intensity  I  ,  and  the  absorption  coefficient  a(s)  all  become  infinite 

9 

for  the  hard  boundary  at  k(s)  *  0,  which  corresponds  to  grazing  reflection 
or  scattering.  This  is  shown  by  case  b  in  Tables  II  and  III.  Furthermore 
ff(s,s')  becomes  infinite  for  the  soft  boundary  at  w  ■  0,  as  is  shown 
in  Table  III,  case  a. 

Multiple  scattering  beyond  double  scattering,  which  has  been 
ignored  in  the  present  analysis,  is  the  mechanism  which  makes  the  physical 
coefficients  finite.  Mathematically  multiple  scattering  corresponds  to 
terms  of  higher  power  in  e,  all  powers  of  which  must  be  included  to  obtain 
finite  results.  Such  terms  are  taken  into  account  in  certain  other  methods 


of  analysis,  such  as  the  smoothing  method  and  the  equivalent  method 
of  partial  summation  of  the  power  series  in  e . 

We  shall  describe  another  method  for  eliminating  the  singularities 
and  obtaining  finite  results  which  are  physically  meaningful.  It  is 
based  upon  representing  the  reflection  coefficient  C(s)  in  the  following 
well-known  admittance- impedance  form: 


C(s) 


+ 


1  +  Q(s) 
1  -  Q(s) 


(42) 


We  choose  the  sign  as  in  Table  I.  To  determine  the  function  Q(s)  we 

require  that  the  power  series  expansion  of  C(s)  given  by  (42)  agree 

3 

with  our  previous  result  (19)  up  to  0(e  ).  This  yields 


Q(s)  -  eh0m(s,s)+e2  J  R(s-s,)ra(s,s,)m(s’ ,s)ds'  +  0(e3)  .  (43) 

The  value  of  Q(s)  for  each  of  the  four  cases  is  shown  in  Table  IV. 

We  shall  call  (42)  the  renormalized  form  of  C(s).  It  does  not 

have  the  defects  of  the  unrenormalized  form  (19).  The  appearance  of  the 

2 

denominator  |l-Q(s)|  in  (41)  when  (42)  is  used  in  it,  suggests  that  we 
should  renormalize  a(s,s')  by  writing 

\  _  4e2R(s-s')|m(s.s,)|2k2(s) 

a(s,s  )  - 5 -  -y — - -  •  fu\ 

1 1-Q(s  * ) | 2  (44) 


This  form  of  a  is  also  free  of  the  defects  of  the  unrenormalized  form  (31) 
When  the  renormalized  C  and  a  are  used  in  the  expression  (41),  it  yields 


the  renormalized  absorption  coefficient  a. 

The  renormalized  expression  (42)  for  C  suggests  that  the  coherent 
field  satisfies  an  admittance  or  an  impedance  boundary  condition  on  the 
mean  surface.  To  show  this,  we  define  the  effective  specific  admittance  of 
the  surface  to  be  Q±1(s)ck(s)/“,  and  we  denote  its  inverse  Fourier  transform  by 
n(x,y,t).  The  sign  is  as  in  Table  I.  Then  we  consider  the  boundary 
condition 

c3  <i)rhp>  +  (2it)  3n(x.y»t)  *  3 *.<4»‘+d»  “0  on  z  -  0.  (45) 

z  t 

Here  *  denotes  convolution  in  x,  y  and  t.  This  boundary  condition 
yields  exactly  the  plane  wave  reflection  coefficient  C  given  by  (42) . 

Alternatively  we  can  introduce  the  reciprocal  of  the  admittance, 
which  is  just  the  effective  specific  impedance  of  the  surface, and  is  given 
by  Q^(s)'o/ck(s)  .  We  denote  its  inverse  transform  M(x,y,t).  Then 
instead  of  (45)  we  write 

+  c(2w)  ^M(x,y,t)  *  dz<^H-<p>  «  0  on  z  *  0.  (46) 

Again  (46)  leads  to  (42).  The  value  of  (45)  and  (46)  is  that  they  can 
be  used  for  a  rough  surface  which  is  not  nearly  flat,  provided  that  its 
radius  of  curvature  is  large  compared  to  the  correlation  length  of  the 
roughness  and  to  the  wavelength. 

There  are  plane  wave  solutions  of  (1)  which  satisfy  (45)  or  (46) 
with  no  incident  wave,  and  which  decay  with  increasing  distance  z  away 
from  the  surface.  Substitution  of  a  plane  wave  into  (45)  or  (46) 
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yields  the  condition  Q(s)  *  1.  This  shows  that  such  waves  correspond 

to  poles  of  the  renormalized  reflection  coefficient  C(s)  given  by  (42). 

The  roots  of  Q(s)  -  1  have  been  examined  for  the  hard  surface,  with  Q(s) 

given  by  case  b  in  Table  IV,  by  Wenzel  (3  ),  and  for  this  and  other 

cases  by  Watson,  Stickler  and  Keller  (4  ) . 

The  renormalization  method  is  arbitrary  in  the  sense  that  the  form 

(42)  of  C(s)  is  postulated.  Once  this  is  done  however,  the  value  of  Q(s)  is 

2 

determined  through  terms  of  order  e  .  Similarly,  the  form  (44)  of 

o  is  specified.  As  a  consequence  the  renormalized  absorption  coefficient 

for  motionless  hard  and  soft  surfaces  vanishes  through  terms  of  order 
2 

e  ,  instead  of  vanishing  identically. 


7.  Embossed  surfaces  and  Twer sky’s  theory 

Some  rough  surfaces  can  be  viewed  as  flat  surfaces  with  random  collections 
of  bumps  or  bosses  glued  onto  them.  This  is  the  model  of  a  rough  surface 
employed  by  Twersky  [  1  ]  in  his  extensive  investigations.  Because  it  is 

so  different  from  that  employed  in  the  application  of  regular  perturbation 
theory,  the  relation  between  perturbation  theory  and  Twersky' s  theory  has 
never  been  adequately  explained.  We  shall  explain  it  by  applying  our  results 
to  embossed  planes  and  then  comparing  them  with  Twersky' s  results.  In  doing 
so  we  shall  extend  the  model  to  moving  surfaces  by  permitting  the  bosses 


to  move. 


We  begin  by  writing  the  equation  of  the  surface  as  z*eh(x,y,t)  with 
N 

h(x,y,t)  -  ^  b[(x-xi)cos0  -  (y-y^sin©^  (y-y^cosS  +  (x-x^sinS^  .  (47) 


Here  b(x,y)  determines  the  shape  of  a  single  boss  located  at  the  origin, 

and  N  is  the  total  lumber  of  bosses. 


All  of  the  bosses  are  of  the  same  shape,  the  i  one  having  been  rotated 


through  the  random  angle  0^  and  translated  to  the  random  position  r^  *  (x^,; 
at  time  t.  Thus  the  statistics  of  the  surface  are  specified  by  the 

joint  probability  distribution  of  all  the  r^  ■  (x^.y^)  and  9^,  which 
are  random  functions  of  t.  However  for  our  purposes  it  will  suffice  to 
specify  the  probability  distribution  for  one  or  two  bosses  at  two  times, 
which  we  shall  now  do. 

Let  us  denote  by  rj,9j,t')  the  probability  density  for 

boss  i  to  be  at  r^  with  orientation  0^  at  time  t,  and  for  boss  j  *  i 
to  be  at  r^*  with  orientation  0^  at  time  t'.  Similarly  let  p^r^.S^t; 
r^,0^,t')  be  the  probability  density  for  boss  i  to  be  at  r  with 
orientation  0^  at  time  t,  and  to  be  at  r^  with  orientation  0^  at 
time  t'.  For  distributions  statistically  stationary  in  space  and  time,  these 
densities  must  have  the  forms 

p2  -  A_2w2(r1-rj,t-t,,0i,0p  r±  and  r^  in  D,  (48) 

•  A-1w1(ri-r^  ,t-t',0i,0p  r±  and  rj^  in  D  .  (49) 

Here  D  is  the  domain  throughout  which  the  bosses  are  distributed  and  A 
is  its  area.  Ultimately  we  will  let  D  become  the  whole  plane  with  A 
and  N  tending  to  infinity  with  a  fixed  ratio  N/A  ■  v.  Thus  v  is  the 

average  number  of  bosses  per  unit  area. 

To  calculate  <h>  =  hQ  we  need  the  probability  density  for  particle  i  to 
be  at  r^  with  orientation  9^  at  time  t,  which  we  can  obtain  by  setting 
t*t*  in  p^.  From  the  definition  of  p^  it  follows  that  at  t*t',  w^  in 
(49)  must  have  the  form  w^r^r^, 0,9^9^)  ■  6(ri-rp6(@i-0pf (9^)  for  r ^ 
in  D.  Then  from  (47)  we  get 


hQ  -  <h>  -  NA_xJjjbidx1dyif(ei)d0i 
D 

Here  b^(r)  is  the  summand  in  (47),  and 
to  x  and  y,  is  the  volume  of  a  boss. 
r(p),  the  auto- cor relation  function  of 

r(p)  -  <h(r+r',t+t')h(r',t,)>  -  hg  = 


w  I  f(8^)d0^  =  w  .  (50) 


v,  the  integral  of  b  with  respect 
Next  we  use  (47 )— (49)  to  calculate 
h,  as  follows: 

2rr2rr 


J 

0  0  DD 


[v  w2(ri-r^,t,0i,0^) 


(51) 


2  2 


+  w1(r1-rj,t,ei,0j)]b1(r+r,)bj(r,)dridrjd01d0^.-v  v 


From  (51)  we  can  calculate  the  spectral  power  density  R(s),  which  is 
the  Fourier  transform  of  r(p),  by  noting  that  the  integrals  in  (51)  are 
convolutions  with  respect  to  r^  and  r  ^ .  Then  we  obtain 


2tt2t r 

R(s)  -  (2t r)4  |  |  tv2W2(s,0i,8j)+vW1(s,8i,0j)]B1(a, 0)^(0, e)d0id0j-v2v26(s),  I 
0  0 

where  B^a, 8)  ■  B[acos8^  +  SsinO^,  -  asin0^  +  ScosS^] .  Here  B  is  the  two 
dimensional  transform  of  b(x,y).  The  results  (50)  for  h^  and  (52)  for 
R(s)  can  be  used  in  (42)-(44)  to  get  C  and  a. 

The  Integrals  in  (52)  simplify  when  the  bosses  are  similarly  aligned, 
when  the  bosses  have  uniformly  distributed  orientations,  and  when  they  are 
axially  symmetric.  For  example  in  the  first  case,  when  all  the  bosses  are 
permanently  aligned  in  the  direction  0*0,  we  take  * 

Wjt(s)5(0i)5(0j),  k  -  1,2.  Then  (52)  becomes 

R(s)  -  (2ir)4[v2W2(s)+vW1(s)]|B(a,3)|2  -  v2v26(s)  .  (53) 


We  shall  now  compare  our  results  for  an  embossed  surface  with  an  old 


result  of  Twersky  [  1].  His  results  apply  to  permanently  aligned  bosses, 


so  we  must  use  (53)  for  R(s).  He  also  considers  small  number  density  v, 

2 


so  we  must  omit  the  v  term  from  (53).  Furthermore  his  surfaces  are  not 
moving,  so  we  must  set  (2ir)2W^(s)  =  <5(uj).  With  these  simplifications,  (53) 
becomes 


R(s)  -  (2Tr)2v|B(a,8)j25(oj)  +  0(v2)  . 


(54) 


When  we  use  (54)  and  (50)  in  (43)  and  (44)  we  obtain 


Q(s)  *  evvm(s,s)+e2(2ir)2v  [  |B(a-a'  ,8-$')  |2m(s,s,)m(s,,s)6(u>-u),)ds, 

1  o  J  (■ 


+  0(e3)  +  0(v2) 


|x-Q(s*)|z 


In  Appendix  B  we  calculate  to  second  order  in  e  the  scattering 
amplitude  f(u,e)  of  a  single  boss  on  a  plane,  for  the  incident  wave 
vector  e  and  the  scattered  wave  vector  u  .  We  shall  denote  the  result 


fg  because  it  is  the  second  Born  approximation.  It  is  given  by  (B6) .  In 


terms  of  f  ,  we  can  write  (55)  and  (56)  in  the  forms 

D 


Q(3) 


wfg(8,e) 


(w/c)k(s) 


(57 


OKs.s’) 


(0)/c)‘ 


fAu>e) 


l-Q(s') 


5  (co-to1)  . 


(58 


Here  §  ■  (e^, e,;, ,-<;,) ,  while  the  relation  between  s  and  s',  and  u  and 
e  is  given  by  (B7)  and  (B8) . 


These  results  (58)  and  (59)  are  of  exactly  the  same  form  as  those 

in  the  abstract  of  Twersky  [1]  with  fg  instead  of  the  exact 

scattering  amplitude  f.  Therefore  the  present  results  for  the  renormalized 

reflection  and  differential  scattering  coefficients  reduce  to  those  of 

* 

Twersky  when  they  are  specialized  to  hard  and  soft  embossed  planes  in 

which  the  bosses  are  permanently  aligned,  not  moving,  and  of  low  number 
density.  Then  they  yield  Twersky's  results  with  f  replaced  by  fg. 

When  f  is  not  a  good  approximation  to  f,  as  is  the  case  for  bosses 
high  compared  to  a  wavelength  or  bosses  having  steep  slopes,  Twersky's 
results  are  more  accurate  than  the  present  ones.  However  the  present 
results  include  surfaces  which  are  not  embossed  planes,  which  are  moving, 
which  may  have  impedance  or  admittance  boundary  conditions,  and  which  can 
have  a  high  density  of  bosses.  In  addition  they  yield  the  two-point  two- 
time  correlation  function  of  the  field*  not  just  the  coherent  field  and  the 
intensity. 

Twersky  [ 5  ]  has  extended  his  theory  to  non-aligned  bosses  and  to 
surfaces  with  a  high  density  of  bosses.  We  plan  to  compare  the  present 
results  with  these  newer  results  of  his. 
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Appendix  A.  Solution  for  A_(s) . 


To  find  As(s)  for  the  case  of  a  soft  boundary,  we  take  the  Fourier 


transform  of  (2a)  to  get 


ico[<f(s,0)  +  f(s,0)]  +  ie  H(s-s')&' [3  $(s',0)  +  3  ’f(s',0)]ds' 


2 

ij-  H(s-s')  H(s'-s")w"[32<i>(s",0)  +  3^(s",0)]ds"ds' 


♦  0(e)  =  0  . 


Then  we  use  (8)  and  (9)  for  Ag  and  Ai  in  (Al)  and  obtain 


ioi[A  (s)  +  A.  (s) ]  +  e  [  H(s-s')a)'k(s')  [A  (s’)  -  A.(s')]ds' 

SI  I  SI. 

-  H(s-s')  |  H(s'-s")co"k2(s") [As(s")  +  Ai(s")]ds"ds* 


+  0(e  )  =  0  . 


Now  we  expand  in  powers  of  e  in  the  form 


A  (s)  =  A(0)(s)  +  eA(1)(s)  +  e2At2)(s)  +  0(e3)  . 


Upon  using  (A3)  in  (A2),  and  equating  to  zero  the  coefficient  of 


each  power  of  e,  we  get  the  equations 


ico[A^(s)  +  A.(s)]  =  0  , 


iwA^(s)  +  H(s-s,)co,k(s')  [A^°')  (s')  -  A.(s')]ds'  =  0  , 


icoA^2)(s)  +  H(s-s,)a)'k(s,)A^1^(s,)ds’ 


-  j  H(s-s’)  H(s'-s">"k2(s")[A^0)(s")  +  Ai(s")]ds»ds"  =  0 


From  (A4)  we  get 


(si  =  -A.  (s)  . 
s  i 


(A7) 


Then  (A5)  yields 


*>> 


=  2ioj 


-1 


HCs-s’lw’kCs'lA^Cs'lds' 


(A8) 


Next  (A6)  has  the  solution 


A<2)(.)  =  -2o)-1 

s 


H(s-s')k(s') 


H(s'-s")o)"k(s")Ai(s")ds"ds'  .  (A9) 


By  substituting  (A7)-(A9)  into  (A3)  we  obtain  the  expansion  of 
As(s).  This  solution  is  of  the  form  (11)  with  the  minus  sign, 
and  with  m(s,s')  =  (iai'/oOkCs1) . 


Appendix  B„  Scattering  Amplitude 

We  consider  a  boss  located  on  the  xy-plane  and  insonified  by  an  incident 
plane  wave  in  the  direction  of  the  unit  vector  e  =  [e^je^.e^],  where  e^  <  0. 
The  incident  field  is  given  by 


iKx,y,z,t)  ■  exp[(27Ti/X)  (e1x+e2y+e3z-ct)] ,  (Bl) 

where  X  is  the  wavelength.  The  scattered  field  <£  has  the  following  expansion 
for  field  points  far  from  the  boss 

<f>(x,y,z,t)  -  ±exp[  (2iri/X)  (e^x+e^y-e^z-ct)  ]  ±  f(u,e)  h0(2TT[x2+y2+z2]i5/X)  • 
exp(-27tict/X)  +  0(X/[x2+y2+z2]!j)  (B2) 


The  plus  sign  applies  to  a  hard  boundary  and  the  minus  sign  to  a  soft  one. 


The  first  term  in  (B2)  is  the  specular  reflection  of  an  incident  plane  wave 
from  the  xy-plane.  The  second  term,  which  is  the  wave  scattered  from  the 
boss,  is  proportional  to 


hQ(r)  =  exp(ir)/ (ir)  ,  (B3) 

2  2  2  W 

and  to  the  scattering  amplitude  f(u,e).  Here  u  =  [x,y,z]/[x  +y  +z  ]2  is 
the  unit  vector  in  the  direction  of  observation. 

If  the  height  of  the  boss  is  small  compared  to  a  wavelength,  and  its 
slope  is  small,  then  a  Born  expansion  can  be  used  to  compute  its  scattering 
amplitude.  We  denote  the  second-order  Born  approximation  of  f  by  f  .  To 

D 

compute  it  we  use  the  formulas  given  in  Sec.  2  with  the  incident  field 
given  by  (Bl)  and  the  surface  shape  h(x,y,t)  =  b(x,y).  From  the  transforms 
of  these  quantitites  and  equations  (7)  and  (9),  we  have 


Ai(s)  =  6(a-2TTe1/X)»6(8-2TTe2/X)6(ort-2TTc/X) 


(B4) 


and 


H(s)  =  B(a,8)6(w)  . 


(B5) 


We  substitute  (B4)  and  (B5)  into  (11)  to  get  an  expression  for  the  amplitude 

of  the  scattered  field  transform  A  (s).  Then  we  obtain  the  transform  of  the 

s 


2.  2.  2, 


scattered  field  by  using  this  result  in  (8). 

We  invert  this  transform  for  large  values  of  ( [xN-yN-z*-]'2/!)  by  using 
the  method  of  stationary  phase.  This  yields  an  expansion  of  <j>  of  the  form 
(B2).  Upon  comparing  terms,  we  get 


fg(u,e)  *  4ir(Gd/c)k(s),{eB(a-a',6-8’)m(s,s’)+e 


B  (a-a" ,  6-6”)  B  (a"-oc  ’ ,  6"-6  ’ ) 
m(s, s")m(s", s' )da"d8"}  .  (B6) 


The  right  side  of  (B6)  is  evaluated  with 


W  ■>  -«  ~-»Vv  V  v 


s  -  [a,S,0)j  *  (2tt/X)  [u1,u2,-c]  (B7) 

and 

s’  =  [a'.B'.w']  =  (27rA)[ei,e2,-c]  .  (B8) 

In  (B6),  k(s)  is  defined  by  (7)  and  m  is  given  in  Table  I  for  both  the 
hard  and  soft  boundary  conditions. 


REFERENCES 

1.  Twersky,  V.,  "On  scattering  and  reflection  of  sound  by  rough  surfaces," 

J„  Acoust.  Soc.  Am.  29 _,  209-225  (1957). 

2.  DeSanto,  J.  A.,  "Theoretical  methods  in  ocean  acoustics,"  in  Topics  in 

Current  Physics,  Vol.  8:  Ocean  Acoustics,  edited  by  J.  A.  DeSanto 
(Springer,  New  York  1979),  Chap.  2,  pp.7-77. 

3.  Wenzel,  A.  R.,  "Smoothed  boundary  conditions  for  randomly  rough  surfaces 

J.  Math.  Phys.  15,  317-323  (1974). 

4.  Watson,  J.  G.,  Stickler,  D.  and  Keller,  J.  B.,  "Reflection,  transmission 

and  scattering  by  rough  interfaces,"  1983,  unpublished. 

5.  Twersky,  V.,  "Reflection  and  scattering  of  sound  by  correlated  rough 

surfaces,  "J.  Acoust.  Soc.  Am.  73,  85-94  (1983). 


(a)  fb)  !  (c)  (d) 


TABLE  I.  The  sign  and  the  function  m^s’). 

O)  Soft:  “1- 2ichQkCs)  +  2[eh0k(s)]2  + 2e2k(s)  |  R(s-s')k(s')ds' 

00  Hard:  1+  2ieh0k(s)  -  2[eh0k(s)] 2-2e2k_1(s)  [rCs-s'^Cs') 

* 

[(w/c)2-aa,-SB’]  •  [(a)’  /c)  2-aa' -SB*  ]ds' 

2 

(c)  Admittance:  l+2ehnu)/ck(s)+2[ehna)/ckCs)]  2+ 4^  !£MSlSU.  ds ' 

U  c  k(s)  k(s') 

2  2 

(d)  Impedance:  -l-2eh0ck(s)/a)-2[eh0ck(s)A>] 2-  —  -a)k^s->  d 

TABLE  II.  Reflection  coefficient  C(s)  through  order  e2. 


M 


(a)  Soft:  4e  R(s-s')  (— )  |k(s')|  k  (s) 


(b)  Hard:  4e2R(s-s» )  (u'/c)2  -  eta1  -  B3’]2sgn  k2(s) 


(c)  Admittance:  4e2R(s-s') (w'/c) 2sgn  k2(s) 


(d)  Impedance:  4e2R(s-s')  2 |k(s')  1 2k2(s) 


TABLE  III.  Differential  scattering  coefficient  o(s,s'), 


(a)  Soft:  ieh0k(s)-e2k(s)  RCs-s'JkCs'Jds* 


Hard:  iehQk  (s)  -  e  2k“ 1  (s )  |r(s  -  s ' )  k' 1  (s  ’ )  [  (w/c)  2-aa » -  SB  T  [  (a)  */c)  2-aa  ’  -  36 '  ]  ds 


2 

(c)  Admittance:  ehQw/ck(s)  +  '~2" j  ~  ^k'Cs^  ^s' 


ehoCk(s)/w  +  |  ds. 


TABLE  IV.  Q(s) ,  which  occurs  in  the  renormalized  reflection  and  scattering 


coefficients,  through  order  e  . 


